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Final  Progress  Report  (July  1,  1976  - September  30,  1977) 

Contract  No.  F4 9620-76-C-0001 
Purchase  Request  No.  7T-0031 

PROTECTION  OF  THE  CARDIOPULMONARY  SYSTEMS  AGAINST  THE  INJURIOUS 
EFFECTS  OF  ACCELERATION 

Evaluation  of  effects  of  time  history  of  gravitational- 
inertial  force  environment  on  regional  dynamic  distri- 
bution of  pulmonary  ventilation,  strains,  edema  and 
blood  volume. 

Introduction 

Under  dyneunic  conditions,  the  mechanical  properties  of  the  lung, 
rib  cage,  diaphragm,  and  abdomen  and  the  coupling  between  these 
various  components  of  the  respiratory  system,  play  an  important 
role  in  determining  the  intrapulmonary  distribution  of  inspired 
gas.  It  is  known  that  the  intrapulmonary  distribution  of  inspired 
gas  is  altered  by  such  factors  as  1)  body  position,  2)  lung 
volume,  3)  chest  wall  geometry  and  compliance,  and  4)  increased 
gravitational-inertial  force  environments.  In  addition,  slippage 
of  lobes  along  interlobar  fissures  might  aid  the  lung  in  conforming 
to  differing  thoraco-abdominal  configurations  without  compromising 
regional  ventilation,  particularly  during  postural  changes  and  ex- 
posure to  increased  gravitational-inertial  force  environments. 


Because  of  the  very  fragile  anatomic  structure  of  the  lungs  and 
the  large  intrathoracic  pressure  differentials  generated  by  changes 
in  the  force  environment  due  to  the  extreme  differences  in  specific 
gravity  of  the  gaseous  contents  of  the  alveoli  and  the  parenchyma 
and  blood  surrounding  these  alveoli,  the  lungs  are  the  most  sus- 
ceptible organ  system  in  the  body  to  malfunction  and  structural 
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injury  due  to  acceleration.  Furthermore,  it  is  quite  possible 
that  the  regional  distribution  of  mechanical  stress  within  the 
lung  is  a predisposing  factor  in  the  localization  of  pulmonary 
pathology.  Therefore,  a quantitative  description  of  the  dynamic 
spatial  distribution  of  lung  parenchymal  stresses  and  strains  and 
the  associated  changes  in  dyneimic  shape  and  dimensions  (i.e.,  the 
dyneunic  geometry)  of  the  intact  lungs  and  thorax  is  required  be- 
fore further  progress  can  be  made  in  describing  the  three-dimen- 
sional mechanical  behavior  of  the  intact  lung.  However,  eluci- 
dation of  the  intrinsic  elastic  behavior  of  the  lung  and  its 
relationship  to  the  chest  wall  and  diaphragm  require  studies  of 
the  intact  lung  with  techniques  possessing  sufficient  temporal 
and  spatial  resolution  to  permit  quantitative  measurements  of  the 
dyneunic  changes  in  the  three-dimensional  parenchymal  stresses  and 
strains,  regional  lung  volumes,  and  thoracic  shape  and  dimensions 
along  with  simultaneous  measurements  of  transpulmonary  and  intra- 
vascular pressures. 

This  final  progress  report  for  Contract  No.  F49620-76-C-0001 , 
covering  the  period  July  1,  1976  to  September  30,  1977,  consists 
of  three  parts.  Parts  I and  II  contain  summaries  of  individual 
studies  performed  during  the  period  stated  above  and  Part  III  con- 
sists of  a list  of  publications  from  our  laboratory  which  were 
based  on  work  supported  by  funds  from  this  contract. 

Part  I:  Regional  Parenchymal  Expansion  and  Lung  Mechanics  in  the 
Intact  Canine  Thorax  from  Measurements  of  the  Displace- 
ments of  Roentgenopaque  Parenchymal  Markers  in  Anesthetized 
Dogs 

Experimental  Preparation.  Regional  lung  parenchymal  expansion  in 
the  intact  thorax  can  be  determined  by  recording  the  displacements 
of  radiopaque  metallic  markers  by  means  of  ccxnputer-based  biplane 
roentgen-videographic  recording  system  which  provides  high  temporal 
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(30/sec)  and  spatial  (+  1.0  mm)  resolution  measurements  of  the 
"tagged"  lungs  (2).  Metallic  markers,  which  consist  of  1-mm-dia- 
meter  spheres  (average  weight:  0.01  gm/marker)  and  2-mm-long 
(1-mm-diameter)  rods  (average  weight:  0.02  gm/marker),  were 
implanted  percutaneously  under  fluoroscopic  control  throughout  the 
parenchyma  of  the  right  lung  of  anesthetized  dogs.  A three-dimen- 
sional grid-like  array  of  30-40  markers  (Figure  1)  is  obtained 
without  causing  physiologically  significant  pneumo-  or  hemothorax 
(2,8).  A minimal  recovery  period  of  2 to  3 weeks  is  provided 
before  the  dogs  are  used  in  an  experiment.  Then,  by  means  of  the 
biplane  videoroentgenographic  system,  it  is  possible  to  continuously 
track  the  displacements,  in  three-dimensional  space,  of  these 
radiopaque  parenchymal  markers  throughout  a complete  respiratory 
cycle  or  maneuver  (e.g.,  a static  pressure-volume  curve). 

DORSAL-VENTRAL  AND  LATERAL  X-RAYS  OF  DOG'S 
THORAX  WITH  LUNG  PARENCHYMAL  MARKERS 


Dortal-Vwtral  Lat»ral 


Figure  1 Dorsal-ventral  and  lateral  x-rays  of  dog's 
thorax  with  metallic  markers  percutaneously 
implanted  in  a three-dimensional  grid-like 
array  throughout  the  parenchyma  of  the  right 
lung. 

Data  Acquisition.  A schematic  diagram  of  the  biplane  roentgen- 
videographic  system  and  recording  ecpiipment  used  to  quantitate  the 
dynamic  displacements  of  lung  parenchymal  markers  during  various 
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respiratory  maneuvers  is  shown  in  Figure  2.  Each  roentgen-video 
imaging  chain  consists  of  an  x-ray  source,  a 9-inch-diameter 
fluoroscopic  screen  and  image-intensif ier , and  a television 
(Plumbicon)  camera.  The  central  axes  of  the  two  roentgen-video 
imaging  systems  can  be  selectively  aligned  from  approximately 
25-30®  (i.e.,  the  angle  of  separation  of  the  imaging  chain  central 
axes)  for  stereo  recording  to  90®  for  orthogonal  recording.  The 
central  axes  of  the  two  imaging  systems  are  aligned  so  as  to  lie 
in  the  same  transverse  plane  relative  to  the  cephalocaudal  axis 
of  the  animal  and  the  animal  is  positioned  so  that  the  lungs  and 
thorax  are  centered  at  the  intersection  of  the  axes  of  the  roentgen- 
video  imaging  chains. 


DYNAMIC  REGIONAL  LUNG  MECHANICS 
Biplane  Video  Roentgenographic  Analysis 


Figure  2 Diagram  of  biplane  roentgen-videographic  system 
for  generation  and  recording  of  data  for  dynamic 
measurements  of  regional  lung  strains  and  volumes. 
Axes  of  the  two  x-ray-video  recording  systems 
can  be  selectively  aligned  for  stereo  (as  shown) 
or  orthogonal  recordings.  Roentgen-video  images 
and  analog  variables  (pressures,  flows,  ECG)  are 
simultaneously  recorded  on  videotape  or  disc. 
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The  fluoroscopic  image  which  is  projected  onto  the  target  of  the 
television  camera  is  recorded  as  a sequence  of  approximately  250 
horizontal  lines  (1  video  field;  l/60th  sec)  which  is  interlaced  with 
a second  field  of  as  many  lines  to  create  a frame  every  l/30th  of 
a second  containing  approximately  500  lines  per  image.  Then  by 
means  of  special  switching  circuitry  (10) , the  two  video  fields 
are  positioned  in  the  left  and  right  halves  of  the  video  picture 
so  that  biplane  videographic  images  are  displayed  and  recorded 
simultaneously  in  the  same  video  field.  The  sequence  of  30/sec 
biplane  images  are  recorded  on  videotape  (Ampex  VR  1100)  in  absolute 
temporal  synchrony  with  up  to  16  channels  of  analog  information, 
such  as  vascular,  pleural  and  airway  pressures,  velocity  and  volume 
of  ventilation  and  the  electrocardiogram  (ECG)  at  sampling  rates 
up  to  approximately  1000  samples/ sec  per  analog  channel  by  means 
of  a video  amplitude  modulator-demodulator  system  (9) . 

At  the  end  of  an  experiment,  pressure,  volume  and  flow  calibrations 
and  the  reference  (x,y,z)  axes  used  for  spatial  localization  of 
the  markers  are  recorded  on  videotape.  The  reference  axes  are 
obtained  by  placing  precision-milled  1-cm  metal  grids  on  the  face 
of  each  image  intensifier  such  that  the  centers  of  the  two  grids 
form  a locus  with  spatial  coordinates  (0,0,0)  to  serve  as  the  center 
of  the  cartesian  reference  system. 

The  geometry  of  the  roentgen-video  imaging  systems  can  be  arranged 
to  accommodate  1)  anesthetized  dogs  supported  either  head -up  or 
head-down  in  a vertically-oriented  radiolucent  water-immersion 
respirator  (2),  2)  awake  dogs  standing  in  a Pavlovian-type  sling 

positioned  on  the  fluoroscopic  table  top  for  studies  in  spontan- 
eously breathing  awake  animals  (1),  and  3)  awake  or  anesthetized 
dogs  breathing  spontanesouly  in  a horizontally-oriented,  radiolu- 
cent, air-filled,  plethysmograph  (Figure  3)  for  continuous  measure- 
ment of  lung  volumes  and  flows. 
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Data  Processing  and  Analysis.  Analysis  of  the  biplane  video  images 
and  multichannel  analog  recordings  is  accomplished  by  transferring 
the  images  from  videotape  to  a stop-action  video  disc  recorder 
(Ampex  DR-lOA)  for  frame-by-frame  identification  of  parenchymal 
marker  positions.  A specially-designed  computer  station  (Figure  2; 
lower  right  corner)  consisting  primarily  of  1)  an  alphanumeric 
keyboard,  to  communicate  with  the  computer,  2)  an  interactive  video 
cursor  which  is  "mixed"  with  the  video  image  thereby  eliminating 
parallax,  3)  a video  disc  controller  (VRC),  and  4)  a high  per- 
formance video  monitor  (Tektronix  632)  , is  used  by  an  operator  to 
track  all  parenchymal  markers  which  are  individually  identifiable 
on  both  images  of  a single  biplane  video  frame.  Once  the  geometric 
coordinates  of  all  parenchymal  markers  have  been  entered  into 
computer  memory,  the  operator  directs  the  computer  via  the  key- 
board to  transfer  the  coordinates  and  the  analog  data  (pressures, 
flows,  EGG)  recorded  with  that  video  image  from  core  memory  to 
digital  magnetic  tape. 

This  regimen  is  then  repeated  on  a frame-by-frame  basis,  for  any 
selected  time  interval  (i.e.,  time  between  freimes  analyzed),  over 
the  duration  of  the  respiratory  maneuver.  In  this  way,  the  dynamic 
displacements,  in  three-dimensional  space,  of  the  parenchymal 
markers  during  various  respiratory  maneuvers  can  be  determined  with 
sufficient  temporal  (30/sec)  and  spatial  (+1.0  mm)  resolution  to 
permit  quantitation  of  regional  lung  parenchymal  displacements  and 
volumes  in  the  intact  thorax. 


Experimental  Procedure.  Adult  mongrel  dogs,  weighing  8-14  kg  with 
parenchymal  markers  implanted  previously,  were  used  for  these  studies. 
The  dogs  were  anesthetized  with  sodium  pentobarbital  (30  mg/kg  I.V.) 
and  a cuffed  endotracheal  tube  was  inserted  into  the  airway.  Inter- 
mittent positive-pressure  breathing  (end- inspiratory  airway  pressure. 


10-15  cm  H2O)  was  maintained  with  a Bird  respirator  (Bird  Oxygen 
Breathing  Equipment,  Inc.,  Palm  Springs,  CA) . Esophageal  pressure 
(P  ) was  measured  with  a latex  balloon  positioned,  under  fluoroscopic 
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control,  in  the  middle  third  of  the  esophagus.  Transpulmonary 
pressure  (^tp^  determined  by  subtracting  airway  / i.e., 

tracheal,  pressure  from  P (5)  by  means  of  special  electronic 
circuitry.  The  electrocardiogram  (ECG)  was  also  recorded.  An 
indwelling  catheter  was  inserted  into  the  cephalic  vein  of  the 
fore  leg  so  that  additional  anesthetic  could  be  given  intravenously 
without  opening  the  plethysmograph  and  disrupting  the  measurements. 

The  dog  was  then  positioned  in  the  plethysmograph  (Figure  3) , in 
either  the  prone  or  supine  position,  and  all  pressure  and  intra- 
vascular catheters,  along  with  the  ECG  leads,  were  exteriorized 
through  air-tight  ports.  The  plethysmograph  was  closed,  sealed 
and  then  aligned  in  the  biplane  roentgen-video  imaging  system 
(Figure  3).  The  animals  breathed  through  a pneumotachograph  (Fleisch 
#1,  Dynasciences  Medical  Products,  Blue  Bell,  PA)  to  measure  velocity 
of  air  flow  and  volume  (integral  of  flow  velocity)  of  air  moved 
during  various  respiratory  maneuvers.  Airway  (P^^)  and  transpul- 
monary (P^p)  pressures,  velocity  and  voliame  of  air  flow  and  ECG 
are  recorded  on  videotape  and,  with  the  exception  of  the  ECG,  on 
a 4-channel  oscillographic  recorder  (Model  7414A,  Hewlett-Packard  Co., 
Palo  Alto,  CA) . Validyne  differential  transducers,  used  to  record 
airway  (DP-9),  and  esophageal  (DP-9)  pressures  and  the  pressure 
drop  across  the  pneumotachograph  (MP  45-2)  were  each  interfaced  with 
the  oscillographic  recorder  via  carrier  amplifiers  (CD-19,  Validyne 
Engineering  Corporation,  Northridge,  CA) . 

Static  pressure-volume  (P-V)  curves  were  obtained  for  the  intact 
and  excised  lungs  by,  starting  from  functional  residual  capacity 
(FRC)  (P^p  = 2 H2O),  manually  inflating  the  lungs  with  a cali- 
brated 1000-ml-syringe  in  a step-wise  manner  (usually  100  ml/volume 
step)  to  total  lung  capacity  (TLC)  (P^p  " then  de- 

flating in  the  same  incremental  manner  back  to  FRC.  Care  is  taken 
to  impose  the  scune  volume  history  on  the  lungs  immediately  prior 
to  doing  the  P-V  maneuvers. 
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Figure  3 Cylindrical,  air-conditioned,  radiolucent  body 
plethysmograph  positioned  in  the  biplane  video- 
roentgenographic  recording  system,  Plethysmo- 
graph permits  quantitation  of  lung  volumes  and 
can  be  used  with  anesthetized  spontaneously 
breathing  or  ventilated  dogs  in  various  posi- 
tions (prone,  supine,  right  and  left  lateral 
decubitus) . 

A.  Effect  of  Cardiogenic  Motion  on  Intermarker  Distance  Measurements. 
The  parenchymal  markers  are  essentially  randomly  distributed 
throughout  the  right  lung  and  hence,  some  are  in  closer  proximity 
to  the  heart  than  others.  It  is  quite  evident  from  viewing  recordings 
of  different  respiratory  maneuvers,  that  parenchymal  markers 
located  near  the  heart  undergo  additional  (i.e.,  unrelated  to 
changes  in  lung  volume)  oscillation  and  displacement  due  to  cardio- 
genic motion.  This  "cardiogenic"  contribution  to  regional  parenchymal 
marker  displacement  is  observable  at  all  lung  volumes  but  its 
magnitude  increases  as  lung  volume  decreases,  i.e.,  as  the  lung 
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becomes  less  stiff.  From  the  physiologic  standpoint,  studies 
indicate  that  this  cardiogenic  motion  contributes  to  gas  mixing  in 
the  lung  and  thereby  reduces  anatomic  dead  space  (3,4) . 

To  study  the  effect  of  cardiogenic  motion  on  regional  parenchymal 
marker  displacements  and  intermarker  distance  measurements,  an 
anesthetized  prone  dog  was  passively  inflated  in  a step-wise  manner 
from  FRC  to  TLC  and  deflated  back  to  FRC.  Biplane  roentgen-video 
images  of  this  P-V  maneuver  were  recorded  on  videotape,  transferred 
to  a stop-action  video  disc,  and  the  geometric  positions  of  the 
markers  were  determined,  as  described  above,  by  two  trained  operators, 
each  repeating  the  process  three  times. 

Initially,  marker  positions  were  determined  at  each  lung  volume 
without  regard  to  the  phase  of  the  cardiac  cycle  at  which  the 
measurements  were  made.  Figure  4 (upper-half)  shows  the  variability 
(i.e.,  standard  deviation)  of  the  intermarker  distances  as  a func- 
tion of  mean  intermarker  distance.  Data  is  based  on  six  deter- 
minations of  intermarker  distances  between  4 marker-pairs  at  each 
of  15  lung  volumes  (N  = 60) . Note  the  small  standard  deviation 
(SD)  of  the  intermarker  distance  measurements  between  markers  close 
together  (15-25  mm)  in  contrast  to  the  relatively  larger  SD  of  the 
measurements  of  distance  between  markers  further  apart  (35-60  mm) . 

The  reason  for  this  is  that  the  effect  of  cardiac  motion  on  marker 
position  diminishes  with  increasing  distance  of  the  marker  from  the 
heart.  For  marker-pairs  which  are  near  the  heart,  cardiogenic 
motion  appears  to  displace  each  marker  approximately  equally  and 
thus  the  distance  between  this  marker-pair  is  little  influenced  by 
the  beating  heart.  Similarly,  for  marker-pairs  which  are  close 
together,  but  distant  from  the  heart,  intermarker  distances  do  not 
appear  to  be  influenced  by  cardiogenic  motion  due  to  the  limited 
transmission  of  this  motion  through  the  lung  parenchyma.  On  the 
other  hand,  for  markers  further  apart  such  that  one  is  near  the  | 

heart  and  one  is  relatively  distant  from  the  heart,  there  is  a | 

greater  effect  of  cardiac  motion  on  intermarker  distance  measure- 
ments. 
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EFFECT  OF  CARDIOGENIC  MOTION  ON  ACCURACY 
OF  INTERMARKER  DISTANCE  MEASUREMENTS 
(6  Trials,  Z Observers ) 
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Figure  4 Effect  of  cardiogenic  motion  on  intermarker 

distance  measurements.  Top-half:  SD  of  inter- 
marker distance  measurements  as  a function  of 
mean  intermarker  distance.  Data  based  on  6 
determinations  of  intermarker  distances  of  4 marker- 
pairs  at  each  of  15  lung  volumes  (N  = 60)  without 
regard  to  phase  of  cardiac  cycle.  Bottom  half: 
Similar  plot  obtained  when  marker  positions  were 
determined  for  each  volume  at  the  same  point  in 
the  cardiac  cycle. 

Figure  4 (lower-half)  is  a similar  plot  obtained  when  the  positions 
of  the  same  markers  were  determined  in  the  same  manner  as  described 
above,  but  at  identical  points  in  the  cardiac  cycle.  Note  the 
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reduction  in  variability  of  the  intermarker  distance  measurements. 
There  is  no  longer  any  effect  of  intermarker  distances  on  repro- 
ducibility. When  cardiogenic  motion  is  taken  into  account,  the 
largest  coefficient  of  variation  observed  was  < 3%  for  a marker 
pair  separated  by  16.5  mm.  It  is  important  to  note  that  movement 
of  the  animal  per  se  does  not  affect  the  measurements  made  with 
this  technique  since  intermarker  distances  (i.e.,  distance  between 
markers  in  three-dimensional  space)  are  computed  rather  than 
displacements  of  markers  from  some  fixed  coordinate  origin. 

Figure  5 shows  the  frequency  distribution  of  the  maximum  differences 
in  intermarker  distance  measurements  observed  out  of  the  6 deter- 
minations for  each  of  the  4 marker-pairs  at  the  15  volumes  examined. 
The  maximum  difference  in  computed  intermarker  distance  was  less 
than  0.7  mm  in  50%  of  the  observations  and  was  less  than  1.0  mm 
in  80%  of  the  observations.  Since  these  measurements  can  be  made 
as  frequently  as  30  times  per  second,  this  biplane  roentgen-video- 
graphic  technique  is  capable  of  providing  temporal  and  spatial 

resolutions  which  are  many-fold  greater  than  that  obtainable  with 

133 

radioactive  gas  (e.g..  Xenon)  techniques  for  measuring  regional 
lung  volumes. 

The  data  presented  in  Figures  4 and  5 were  obtained  from  a static 
pressure-volume  maneuver,  i.e.,  a step-wise  inflation-deflation 
from  FRC  to  TLC  and  back  to  FRC.  From  such  a procedure,  it  is 
possible  to  plot  a pressure-volume  curve  which  provides  limited 
information  regarding  the  overall  mechanical  characteristics  of 
the  lung  in  the  intact  thorax.  Since  the  biplane  roentgen-video- 
graphic technique  is  capable  of  accurately  determining  the  spatial, 
i.e.,  the  three-dimensional,  coordinates  of  each  parenchymal  marker, 
it  is  possible  to  construct  regional  pressure-intermarker  distance 
curves  for  different  lung  regions.  Figure  6 shows  two  such  curves 
for  marker-pairs  located  in  the  upper  (marker  pair  2,3)  and  lower 
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Figure  5 A plot  of  the  frequency  distribution  of  the 
maximum  differences  in  intermarker  distance 
measurements  obtained  by  2 operators  (3 
trials  each)  using  an  operator-interactive, 
computer-generated  video  cursor  to  identify 
the  spatial  coordinates  of  8 markers  at  each 
of  15  volumes  (N  = 60)  during  a static  pres- 
sure-volume maneuver.  Measurements  at  each 
lung  volume  were  made  at  the  same  point  in  the 
cardiac  cycle. 

(marker  pair  5,7)  regions  of  the  right  upper  lobe  of  the  intact 
canine  lung.  The  data  are  the  means  + SD  of  6 separate  determina- 
tions by  2 operators  (3  trials  each)  for  the  two  marker-pairs  at 
each  of  15  lung  volumes.  The  variability  (i.e.,  SD)  of  these 
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measurements  is  quite  small  as  evidenced  by  the  vertical  (SD)  bars. 
Note  also  the  hysteresis  between  the  inflation  (•)  and  deflation  (0) 
portions  of  the  maneuver  as  is  seen  in  the  conventional  whole  lung 
P-V  curve. 

If  the  lung  expanded  in  a uniform  or  isotropic  manner,  i.e.,  for  a 
given  change  in  stress  (transpulmonary  pressure) , the  fractional 
elongation  of  the  lung  parenchyma  along  the  x,  y,  and  z axes  was 
equivalent,  it  would  be  possible  to  calculate  the  spatial  distri- 
bution of  changes  in  regional  lung  volumes  by  simply  cubing  the 
distances  between  marJcer-pairs  located  in  the  various  regions  of 
the  lung.  However,  it  has  been  shown  that  lung  expansion  is  not 
isotropic  and  that  strains  measured  between  two  markers  will  depend 
on  the  spatial  orientation  of  the  markers  within  the  lung  parenchyma 
(2).  Hence,  pairs  of  markers  can  be  used  to  quantitate  inter-  and 
'ntralobar  strains,  but  cannot  be  used  to  calculate  dynamic  changes 
in  regional  lung  volume  (i.e.,  ventilation).  However,  the  latter 
can  be  quantitated  by  selecting  sets  of  4 non-coplanar  markers, 
with  any  given  set  of  four  contained  within  the  same  lobe,  to 
define  the  vertices  of  a tetrahedron.  Then  by  simply  calculating 
the  volume  of  each  tetrahedron  directly  from  solid  geometry,  the 
spatial  and  temporal  distributions  of  changes  in  regional  lung 
volume,  relative  to  the  lung  itself  (in  contrast  to  inspired  radio- 
active gas  techniques)  can  be  obtained  throughout  a respiratory 
maneuver.  The  dynamic  changes  in  relative  position  of  the  markers, 
and  hence  the  volume  descrilaed  by  them,  can  thf-n  be  further  analyzed 
by  well-established  methods  in  mechanics  (e.g.,  the  affine  trans- 
formation) . Software  to  implement  the  affine  transformation  has  re- 
cently become  operational  and  initial  results  are  currently  being 
analyzed. 
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INTERMARKER  DISTANCE  - PRESSURE 
CURVES  FOR  INTACT  CANINE  LUNG 

(6  Trials,  2 Observers) 


TRANSPULMCNARY  PRESSURE,  cm  HgO 


Figure  6 Intermarker  distance  - transpulmonary  pressure 
curves  obtained  for  the  intact  canine  lung. 
Intermarker  distance  measurements  represent 
the  means  + SD  of  6 separate  determinations  by 
2 operators  (3  trials  each)  during  a static 
pressure-volume  maneuver.  Measurements  at 
each  lung  volume  were  made  at  the  scune  point 
in  the  cardiac  cycle. 

B.  Regional  lung  expansion  in  the  intact  vs.  excised  canine 
lung  at  TLC.  To  analyze  regional  lung  deformation  with  the 
goal  being  to  delineate  the  factors  responsible  for  nonuniform 
ventilation,  it  is  important  to  know  if,  at  total  lung  capacity 
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(TLC),  the  lung,  in  the  intact  thorax,  is  deformed  from  its 
shape  under  conditions  of  uniform  or  isotropic  expansion.  Studies 
of  dog  lungs  quick  frozen  at  high  inflation  pressures  reveal  no 
systematic  regional  differences  in  alveolar  size  (6,7).  This 
suggests  that  at  TLC  the  lung  is  uniformly  expanded. 

To  compare  regional  lung  expansion  in  the  intact  and  excised  lung 
at  TLC,  we  selected  four  dogs  with  good  three-dimensional  distribu- 
tions of  parenchymal  markers.  Static  pressure-volume  curves  were 
obtained  for  the  intact  lungs  by  placing  the  dogs  in  the  plethysmo- 
graph  (Figure  3)  in  the  prone  position  and  inflating  the  lungs  in 
a step-wise  manner  from  FRC  to  TLC  and  then  deflating  in  the  same 
incremental  manner  back  to  FRC.  The  animals  were  then  sacrificed 
and  exsanguinated  and  the  lungs  were  excised  through  a median 
sternotomy.  No  pleural  adhesions  were  encountered  in  any  of  the 
animals  studied.  The  entire  excised  lung  was  then  returned  to  the 
plethysmograph  and  placed  in  the  sling  in  as  close  to  the  same 
geometric  orientation  as  it  assumed  in  the  intact  prone  animal  as 
was  possible.  Static  P-V  curves  were  again  obtained  by  inflating 
the  lung  from  FRC  (P^p  ~ H2O)  to  TLC  ~ 25  cm  H2O)  and 

then  deflating  back  to  P^p^  2.5  cm  H2O  in  a step-wise  manner. 

Biplane  roentgen-video  images  of  these  maneuvers  were  recorded  on 
videotape  and  the  distances  between  all  intralobar  marker-pairs 
were  determined  at  total  lung  capacity,  as  described  above. 

The  results  from  the  four  dogs  are  shown  in  Figures  7 and  8.  The 
distances  between  all  intralobar  marker-pairs  at  TLC  are  plotted 
for  the  intact  prone  dog  (ordinate)  versus  the  excised  lung  (abscissa) 
at  the  same  inflation  pressure.  A good  correlation  (range; 
r = 0.966  to  0.999)  and  a high  degree  of  linearity  of  the  distances 
between  intralobar  marker-pairs  located  in  all  three  lobes  of  right 
lung,  intact  and  excised,  was  observed.  Note  that  the  regression 
line  is  very  close  to  the  line  of  identity.  Since  the  markers 
were  essentially  randomly  distributed  in  the  lung  parenchyma,  this 
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significant  linear  correlation  indicates  that  the  relationships 
of  intralobar  intermarker  distances,  in  the  intact  and  excised 
lung,  was  independent  of  spatial  orientation  and  location  within 
the  lung.  The  close  correspondence  between  intralobar  intermarker 
distances  indicates  that  each  lobe  was  not  only  at  the  scime  volume 
but  was  expanded  in  an  identical  manner  intact  and  excised.  Assuming 
that  at  TLC,  i.e.,  at  high  inflation  pressures,  the  excised  lung  is 
stiff  enough  not  to  be  deformed  by  its  own  weight,  the  lung  should 
be  uniformly  or  isotropically  expanded.  Therefore,  changes  in 
regional  parenchymal  shape  and  dimensions  from  that  pertaining  at 
TLC  represent  changes  from  the  state  of  uniform  expansion. 


INTERMARKER  DISTANCES  AT  TLC  IN 
INTACT  vs.  EXCISED  CANINE  LUNG 

(P,p  » 25  cmHzO  at  TLC) 


INTERMARKER  DISTANCE  at  T LC EXCISED  • 


Figure  7 Comparison  of  distances  between  intralobar 
marker-pairs  in  the  right  lung  of  2 dogs 
at  total  lung  capacity  (TLC)  with  the  lungs 
in  the  intact  thorax  (ordinate)  and  excised 
(abscissa) . 
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INTERMARKER  DISTANCE  AT  TLC  IN 
INTACT  vs.  EXCISED  CANINE  LUNG 

( Pfp  - 25  cm  HgO  at  TLC ) 


INTERMARKER  DISTANCE  at  T L C EXCISED  • 


Figure  8 


I 

Similar  to  Figure  7 for  2 additional  dogs. 


C.  Regional  lung  expansion  at  TLC  in  the  intact  thorax  in  the  prone 
and  supine  positions.  A similar  series  of  experiments  were  per- 
formed in  intact  anesthetized  dogs  in  the  prone  and  supine  body 
positions.  Three  dogs,  with  good  three-dimensional  marker  distri- 
butions, were  used  for  this  study,  one  of  which  (E-123)  was  used 
in  both  studies.  Each  dog  was  anesthetized  and  positioned  in  the 
plethysmograph  in  the  prone  and  supine  positions  and  pressure- 
volume  maneuvers  were  performed,  as  described  above,  in  both 
positions  for  each  animal.  The  results  from  the  three  dogs  are 
shown  in  Figure  9.  The  distances  between  all  intralobar  marker- 
pairs  were  determined  at  TLC  = 25  cm  H2O)  and  compared  when 

the  animals  were  in  the  supine  (ordinate)  versus  prone  (abscissa) 
position.  Again,  there  is  a significant  linear  correlation  (range; 
r = 0.96  to  0.99)  of  distances  between  intralobar  marker-pairs 
in  all  three  lobes  of  the  intact  right  lung  when  compared  at  the 
scime  transpulmonary  pressure  in  the  supine  vs.  prone  positions. 


I 
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EFFECT  OF  BODY  POSITION  ON  REGIONAL  LUNG  EXPANSION  AT  TLC 

(P,p  -25  cm  HgO  at  TLC) 
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Figure  9 Comparison  of  distances  between  intra-  and 
interlobar  marker-pairs  in  the  intact  lungs 
of  3 dogs  at  total  lung  capacity  (TLC)  with 
each  dog  in  the  supine  (ordinate)  and 
prone  (abscissa)  body  positions. 

The  close  correspondence  of  distances  between  intralobar  marker- 
pairs  indicates  again  that  each  lobe  in  the  intact  thorax  was  not 
only  at  the  seune  volume,  but  was  expanded  in  a similar  manner  when 
the  animal  was  in  the  prone  and  supine  positions. 

In  contrast  to  intralobar  intermarker  distances,  distances  between 
interlobar  marker-pairs  (Figure  9),  i.e.,  marker-pairs  that  span 
fissures,  do  not  show  the  same  high  degree  of  linear  correlation 
(range:  r = 0.77  to  0.86)  that  was  observed  with  the  intralobar 
marker-pairs.  These  data  suggest,  therefore,  that  individual  lobes 
shift  their  positions  relative  to  each  other  when  the  position  of 
the  animal  is  changed  from  the  prone  to  the  supine.  Slippage  of 
lobes  along  interlobar  fissures  could  aid  the  lung  in  conforming  to 
differing  thoraco-abdominal  configurations  and  shifts  in  thoracic 
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contents  associated  with  postural  changes  without  compromising 
regional  expansion  of  the  individual  lung  lobes. 


INTERMARKER  DISTANCES  AT  TLC  IN 
ANESTHETIZED  vs.  ANESTHETIZED  - PARALYZED  DOG 
(Nembutal,  Pavulon,  Pfp  ~ 25  cm  HgO  at  TLC) 


INTERMARKER  DISTANCE 
at  T LC/^tiESTHETtZED-PARALYZED  < 


Figure  10  Comparison  of  distances  between  intralobar 
marker-pairs  in  the  intact  lung  at  total 
lung  capacity  in  an  anesthetized  dog  before 
(ordinate)  and  after  (abscissa)  muscle 
paralysis . 

Initial  studies  in  one  dog  indicate  that  regional  lung  expansion 
at  TLC  is  the  same  in  the  anesthetized-paralyzed  as  in  the 
anesthetized  dog  (Figure  10)  which  suggests  that  regional  lung  and 
chest  wall  geometries  are  not  significantly  altered  by  muscle 
paralysis . 
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Thus,  we  have  observed  a highly  significant  linear  correlation  of 
distances  between  intralobar  marker-pairs  in  all  three  lobes  of 
the  right  lung  in  the  intact  thorax  of  anesthetized  dogs  in  the 
prone  and  supine  positions  at  25  cm  H2O  transpulmonary  pressure. 

A similarly  significant  linear  correlation  was  observed  when  the 
lung  in  the  intact  thorax  was  compared  with  the  excised  lung  also 
at  the  same  transpulmonary  pressure  (i.e.,  P^p  = 25  cm  H2O) . These 
data,  therefore,  indicate  that  regional  lung  expansion  at  total 
lung  capacity  is  identical  in  intact  dogs  in  the  supine  and  prone 
positions  and  in  individual  lobes  of  the  excised  lung. 


Generally,  regional  lung  volumes  are  expressed  as  a percentage 
of  regional  volume  at  TLC.  It  is  presumed  that  at  TLC,  all  regions 
of  the  lung  are  expanded  equally.  Morphometric  studies  of  lungs  of 
dogs  quick  frozen  at  high  inflation  pressures  = 30  cm  H2O) 

revealed  no  systematic  regional  differences  in  alveolar  size  (6). 

In  addition,  Hogg  and  Nepszy  (7)  examined  lung  tissue  density  in 
dogs  frozen  head-up  at  TLC  (P^-p  = 30  cm  H2O)  and  found  that  regional 
volumes  are  reasonably  uniform  at  TLC  except  for  the  extreme  upper 
and  lower  lung  regions.  We  are  aware  of  no  previous  studies  of 
regional  lung  expansion  in  the  intact  thorax  of  living  (i.e., 
anesthetized)  dogs  with  direct  comparisons  in  the  same  lungs  after 
removed  from  the  thorax. 


The  results  reported  here  indicate  that  regional  lung  expansion  at 
TLC,  defined  in  these  studies  as  that  lung  volume  obtained  when 
P^p  = 25  cm  H2O,  appears  to  be  identical  in  dogs  with  an  intact 
thorax  in  both  the  supine  and  prone  body  positions  and  in  the 
individual  lobes  of  the  excised  lung.  If  we  ass'-'me  that  excised 
lungs  are  uniformly  expanded  at  TLC  because,  at  high  inflation 
pressure,  and  in  the  presence  of  a uniform  pleural  pressure,  they 
are  too  stiff  to  be  deformed  by  their  own  weight,  then  we  can 
conclude  that  at  TLC,  regional  lung  expansion  is  uniform  or  iso- 
tropic in  the  dog  with  an  intact  thorax.  This  follows  from  the 
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fact  that  since  the  parenchymal  markers  were  randomly  distributed 
throughout  the  lung,  the  close  agreement  in  intralobar  intermarker 
distances,  intact  and  excised,  indicates  that  the  lungs  were  also 
uniformly  expanded  in  the  intact  thorax.  Therefore,  in  analyses 
of  regional  lung  deformation,  changes  in  lung  parenchymal  shape 
and  dimensions  from  those  pertaining  at  TLC  represent  changes  from 
the  state  of  uniform  or  isotropic  expansion.  This  does  not  mean 
that  overall  lung  shape  intact  and  excised  was  the  same  or  even  that 
thoracic  cavity  shape  at  TLC  in  the  supine  vs.  prone  animal  was 
necessarily  the  same.  But  the  data  do  show  ti.at  at  TLC,  while 
overall  lung  shape  may  be  different  in  different  body  positions, 
individual  lobes  remain  uniformly  expanded  and  could  presumably 
accommodate  different  thoracic  shapes  by  altering  their  positions 
relative  to  each  other. 

Part  II:  Finite-Element  Analysis  of  Strain  Variations  in  an 
Excised  Lobe  of  the  Canine  Lung 


Experimental  Preparation.  A 10-kg  dog  was  anesthetized  with  an 
intravenous  injection  of  sodium  pentobarbital  (30  mg/kg)  and  the 
right  lower  lobe  of  the  lung  was  excised  via  a median  sternotomy. 
The  lobe  was  rinsed  with  physiologic  saline  and  kept  moist 
throughout  the  experiment  by  spraying  lightly  with  saline  at 
frequent  intervals.  Metallic  markers,  consisting  of  1-mm  diameter 
spheres  and  cylindrical  rods  of  2-mm  in  length  were  affixed  to  the 
pleural  surface  of  the  excised  lobe.  The  upper-half  panels  of 
Figure  11  shows  two  photographs  of  the  excised  lobe  in  the  position 
which  it  was  studied  and  viewed  from  front  and  rear,  respectively. 

A three-dimensional  drawing  of  the  lobe  depicting  the  positions 
of  the  markers  is  shown  in  the  lower-half  panel  of  Figure  11. 

To  inflate  the  lobe,  a cannula  (6-mm  inside  diameter)  was  in- 
serted into  the  main  bronchus  of  the  lobe  and  tied  securely  in 
place.  Airway  pressure  was  measured  with  a water  manometer  con- 
nected via  a Y-tube  to  the  cannula.  The  transpulmonary  pressure. 
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P^p,  of  the  lobe  was  then  equal  to  the  airway  pressure  since 
pleural  pressure  was  equal  to  ambient  pressure,  except  at  the  base 
where  the  weight  of  the  lobe  was  supported. 


EXCISED  LOBE  OF  CANINE  LUNG  WITH 
INTRAPARENCHYMAL  AND 
SURFACE  METALLIC  MARKERS 

(Dog,  10  kg,  Right  Lower  Lobe) 


3-D  Drawing  of  Lobe 
Showing  Positions 
of  Markers 


Figure  11  Upper-half:  Photographs  of  excised  right  lower 
lobe  of  canine  lung  with  intraparenchyma 1 and 
surface  metallic  markers.  Lower-half:  Three- 
dimensional  drawing  of  lobe  showing  spatial 
locations  of  the  markers. 


Static  pressure-volume  (P-V)  maneuvers  were  performed  with  the 
lobe  in  the  position  shown  in  Figure  11,  by  manually  inflating 
the  lobe  from  P.  = 4 cm  H„0  to  total  lobe  capacity  (TL  C)  de- 
fined  as  P^p  = 20  cm  H2O,  and  then  deflating  back  to  P^p  = 4 cm 
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H2O  in  a stepwise  manner  (50  ml/volume  change)  Biplane  roentgen- 
videographic  images  of  this  P-V  maneuver  were  recorded  on  videotape, 
transferred  to  a stop-action  video  disc  and  the  spatial  (i.e., 
three-dimensional)  coordinates  of  the  markers  were  determined  at 
each  volume  as  described  in  Part  I (above)  and  previously  (2). 

The  three-dimensional  coordinates  of  the  metallic  markers,  re- 
lative to  an  arbitary  cartesian  coordinate  system  are  presented 
in  Table  1 (see  page  26  ) . Four  lobe  volumes  were  studied  beginning 
with  TL  C (P.  = 20  H~0)  followed  by  80%,  60%,  and  40%  of  TL  C at 

which  time  P^^  was  equal  to  4 cm  H2O. 

Formation  of  Tetrahedral  Elements.  A solid  of  complicated  geometry 
to  be  analyzed  by  the  finite-element  technique  has  to  be  divided 
into  a finite  number  of  non-overlapping  elements  of  simple  shape 
(11) . The  most  basic  shape  of  elements  hitherto  used  in  three-dimen- 
sional finite-element  analysis  is  a tetrahedron.  The  partitioning 
of  the  excised  right  lower  lobe  of  canine  lung  into  non-overlapping 
tetrahedrons  can  be  implemented  by  selecting  the  metallic  markers 
affixed  to  the  pleural  surface  of  the  lobe  as  vertices  to  be 
shared  by  adjoining  tetrahedrons.  The  step-by-step  partitioning 
procedures  are  illustrated  in  Figures  12  and  13. 

First,  the  lobe  under  consideration  is  approximated  by  the  model 
displayed  at  the  upper  left  corner  of  Figure  12  consisting  of  a 
space  bounded  by  the  four  lines  1-5-13-18,  2-6-14-19,  2-9-12-20, 
and  3-4-10-16  connecting  the  selected  metallic  markers.  The 
model  is  then  partitioned  into  two  regions,  A and  B,  by  use  of 
the  triangular  cutting  planes  1-2-9,  1-5-9,  5-9-13,  9-12-13, 

12-13-20,  and  13-18-20.  Most  of  the  major  bronchi  are  contained 
in  region  A whereas  region  B is  composed  mostly  of  alveoli.  Such 
partitioning  thus  allows  the  analvsis  of  two  regions  with  rather 
contrasting  anatomic  properties.  To  facilitate  the  analysis  of 
strain  variations  during  deflation  of  the  lobe  in  the  direction 
of  the  vertical  body  axis,  regions  A and  B are  both  further 
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, recorded  by  biplane  videoroentqenographic  technique. 
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divided  into  three  layers,  I (upper),  II  (middle),  and  III  (lower) 
by  use  of  the  triangular  cutting  planes  4-5-9  and  10-12-13  for 
region  A and,  5-6-9  and  12-13-14  for  region  B. 

Three  tetrahedral  elements  are  then  formed  in  each  layer  except 
the  layer  I of  region  B for  which  only  two  tetrahedral  elements 
can  be  formed.  This  results  in  a total  of  seventeen  tetrahedral 
elements,  of  which  the  numbering  and  their  respective  vertices  are 
delineated  in  Figure  13. 

Calculation  of  Normal  and  Shearing  Strains.  For  evaluation  of 
the  strains  in  the  lobe  during  the  time  interval  of  At  = t'-t,  the 
spatial  coordinates  of  the  implanted  markers  recorded  at  the 
instants  t and  t'  can  be  utilized.  The  three-dimensional  dis- 
placement of  any  point  P in  a tetrahedron  formed  with  four 
implanted  markers  as  vertices  can  be  assumed  as  linear  functions  of 
its  coordinates; 


u,,  = 


“l^°‘2''P 


x.,+a  -,yp+a 


4^P 


Vp  = 6i+82Xp+e3yp+64Zp 


Wp  = Yi+Y2Xp+Y3yp+Y4Zp 


(la,b,c) 


Since  the  displacements  of  the  four  vertices  of  the  tetrahedron  must 
also  satisfy  the  above  equations,  it  leads  to  twelve  equations 
which  suffice  to  solve  for  the  twelve  unknown  coefficients  a's,  6's, 
and  y's.  The  results  are,  for  i=l,2,3,4 


“l  “ W Si“k' 


®i  = W ^ki'^k' 
k=l 


ZV  . ^ki'^k 
k*l 


-28- 


where  V is  the  volume  of  the  tetrahedron  ecjual  to  one  sixth  of 
the  determinant  of  the  matrix 


1 x^ 

1 X2  'i2  ^2 

1 X3  y3  Z3 
1 X4  y^  Z4 


(3) 


and  in  Equations  (2)  are  the  cofactors  of  the  matrix  M which 
is  composed  of  the  spatial  coordinates  of  the  four  vertices  of 
the  tetrahedron. 

It  is  well  known  from  the  theory  of  elasticity  that  the  normal 
strains  g and  shearing  strains  y are  related  to  the  displacements 
u,  V,  and  w by  the  equations 


{ e}  = [e^e  e„Y  y y ,T 

X y z ' xy'yz' zx] 


r3u  3v  . ^ + iii  4.  iX 

^9x  9y  9z  9x  9y  9y  9z  9z  9x^ 


In  Equations  (3)  and  (4),  the  notation  { } and  []  are  used  for 
vector  and  matrix,  respectively,  and  T denotes  the  transposition 
of  a matrix.  Since  the  displacements  are  assumed  to  be  linear 
functions  of  the  spatial  coordinates  in  Equation  (la,b,c),  the 
strains  obtained  by  partial  differentiation  of  the  displacements 
are  therefore  independent  of  the  coordinates.  The  strains 
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so  calculated  represent  the  average  values  in  the  entire  tetra- 
hedron. Substitution  of  Equations  (1)  and  (2)  into  Equations  (4) 
gives  the  required  formula  for  the  strain  calculations.  That  is 


where 


{ e}  = IB]{6} 


(5) 


{ 6 } = [u. 


'^l  ^2 


V2  W2 


u. 


^3  '^3  ^^4 


V4  W4] 


(6) 


[B]  = [B^  B2  B3  B^]  (7) 

and  for  i=l, 2,3,4 


^i2 

0 

0 

Ci3 

0 

[B,]  = 

0 

^i3 

0 

Ci2 

Ci4 

0 

0 

Ci4 

0 

Ci3 

Hence,  the  procedure  of  calculating  the  normal  and  shearing  strains 
is  first  to  construct  the  matrix  M and  then  compute  its  cofactors 
j to  form  the  matrix  B,  and  finally  use  Equation  (5)  . Not  only 
are  the  spatial  coordinates  of  the  four  vertices  of  the  tetrahedron 
at  the  last  instant  t'  required,  but  also  the  cha nges  in  their 
spatial  coordinates  are  needed  for  the  calculation  of  strains. 

The  next  section  discusses  how  the  displacements  of  a tetrahedral 
element's  vertices  are  to  be  calculated  based  on  the  changes  in 
their  spatial  coordinates. 

Engineering  Strains  and  Instantaneous  Strains.  Suppose  that  the 
spatial  coordinates  of  the  implanted  markers  are  recorded  at 
instants  t = t^^,  t2,  ...,  t^.  Any  one  of  these  instants  can  be 
chosen  as  the  reference  fraune,  t^.  Both  the  engineering  strains 
and  instantaneous  strains  of  the  tetrahedron  can  be  computed. 
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Engineering  strains  are  determined  using  Equation  (5)  with  the 
displacements  calculated  at  instant  t^  as  follows; 


'i  " 


W.  = Z.  (tj)-Z.  (t^) 


(9a,b,c) 


Equations  (9)  are  to  be  used  for  i=l,2,3,4  to  correspond  to  the 
four  vertices  of  the  tetrahedron. 

Instantaneous  strains  at  instant  t^  are  determined  using 
Equation  (5)  and  the  spatial  coordinates  of  the  vertices  at  two 
consecutive  instants  and  t ^ . The  displacements  are  computed 

with  equations,  for  i=l,2,3,4 


u . 
1 


X. 


Vi  = yi(t3)-yi(tj.i) 


Wi  - Zi(tj)-2i(t._i) 


(10a, b,c) 


Numerical  Results  and  Discussion.  Based  on  the  data  presented 
in  Table  1 of  the  changes  in  the  spatial  coordinates  of  metallic 
surface  markers  recorded  during  deflation  of  an  excised  lobe  of 
canine  lung,  the  engineering  (cumulative)  normal  and  shearing 
strains  calculated  by  the  finite-element  method  are  presented  in 
Tables  2 and  3,  respectively.  These  are  region-by-region, 
tetrahedron-by- tetrahedron  results.  As  mentioned  above,  the 
strain  values  determined  by  the  tetrahedral-element  approach 
represent  the  averages  for  all  of  the  points  contained  in  that 
tetrahedron.  To  improve  the  accuracy  of  the  strain  results, 
their  values  at  the  partitioning  nodes  (selected  metallic  markers) 
can  be  modified  by  averaging  the  respective  strain  values  in 
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Table  2.  Changes  in  engineering  normal  strains  in  the  6 regions  and  17  tetrahedrons  of  the 


32 


I 


03 

to 

‘XU 

x::  — < D3 

■«» 

X— «• 

c:  ro  ft) 

*— < 

*— • 

*5  r-*-v3 

^ -4 

1 cr  - 

--• 

k— 4 

ro  a»  0 

“^3-3  1 

t/>  ro  t/y 

O) 

c 

3 

CTi 

tn 

CO  ro 

0 

CO 

CO 

cn 

cn 

Xu 

to  ro  — * 

— 

i Cl 

1 

(3  O 

1 

O 

o 

1 I 

O CO 

1 1 

0 0 

0 

0 

CD 

w3 

CD 

1 1 

000 

-- 

ro 

ti-i 

to  O 

0 

0 

fO 

ro 

CO 

to 

CO 

C-4  — ' to 

*-4 

o 

o 

tn  •£) 

0 

•0 

OJ 

o-> 

r\j 

cn 

0 

X 

f-j 

CO 

'O 

ay  ro 

C44 

tn 

—I 

••n 

0 

0 

J3 

rj 

— -'CD 

U1 

fv5 

CM 

o 

*•  cn 

‘D 

ID 

fN3 

U1 

:n 

CO  m — * 

1 

I 

1 

1 

1 

, 

1 

1 1 

0 

i 

o 

O 

o 

<3  O 

0 

0 

.3 

0 

•D 

•:d 

0 

0 

000 

y 

UJ 

« « 

cr. 

ro  — • 

rj 

ro 

n 

t 0 

Cu 

0 

0 0 <aJ 

.Ty 

ro 

ay 

x» 

<0 

0 

cn 

ro 

ay 

CO 

<-n  — * to 

*< 

cr»  o 

'O 

— j 

tjy  '-o 

cn 

•D 

cn 

■D 

XU 

vD 

Xu  'vj  ’'O 

tsi 

ro 

Ol 

Ul 

03 

CO  '*0 

CO 

0 

ro 

-X» 

«D 

CO 

to 

cn 

CO  CD  0 

•H 

j 

1 

J 

r- 

, 

, 

. 

, 

, 

, 

, 

, 

0 

1 o o 

O 

0 0 

0 0 

0 0 

0 

000 

0 — • 0 

ro 

Uj 

ro 

0 

rj 

CO 

ro 

(O 

-u 

fn 

0 ro  — * 

-< 

O UJ 

•^4 

CO 

<o 

ro 

ro 

0 

CD 

CD 

0 XU  0 

rx 

cry 

*sj 

CD 

to 

CO 

CO 

CD 

ro 

0 ro  oj 

X 

cn 

ro 

•'O 

-T* 

0 

— j 

ro 

0 

-0 

0 

CD 

*0  X»  0 

■ 

o 

» 

O 

\ 

o 

o 

) t 

0 0 

0 

0 

1 

0 0 

f 

ro 

0 

0 

0 

1 

000 

■ 

o 

CO 

CO 

cn  — • 

0 

n 

— J 

Xu 

CO 

XU 

XU 

— * 0 to 

-< 

iT> 

ro 

0 ^ 

CD 

to 

0 

CO 

CD 

cn 

cn 

'-O  0 

X 

—4 

cn  oy 

CD 

CO 

CO 

cn 

to 

cn 

XU  00  — * 

>< 

— * 

—1 

kO  — • 

CO 

XU 

.D 

CO 

CO 

— * 

ro 

ro 

— ' rrs 

0 

0 

J 

J 

. 

1 <3  0 0 

O <3  O 

0 

0 

0 

CD 

000 

000 

C7> 

cn 

^ cn  o 

cn 

CD 

iu 

0 

XU 

0 

— * ro  — - 

-< 

03 

ro 

a> 

— < cn 

0 

CD 

0 

u> 

0 

0 to 

0 

— * Xu  — • 

*< 

ro 

CO 

0 (O 

ro 

CD 

'O 

fo 

0 

0 

CO 

tn 

ro  XU  — * 

fSI 

C>J 

ay 

ay 

cn  CO 

oy 

CO 

ro 

-TU 

to 

(T> 

00  U3  — » 

-H 

' 

r- 

o o o 

0 0 

0 0 

000 

000 

000 

(_o 

ay 

CO 

r\3  cn 

cn 

IN3 

CO 

•4-» 

CD 

tn 

XU 

cn 

0 ^ CO 

-< 

cjy 

-o 

-fi*  O CO 

'vj 

XU 

cn 

ro 

'O 

to 

03 

03  XU  LO 

ro 

ro 

Oy 

-t*  tjy 

0 

cn 

0 

0 

vD 

CO 

—* 

tn  XU  XU 

X 

— * 

LO 

CO 

ro  CO 

CO 

cn 

cn 

cn 

to 

CO 

cn 

0 

cn  C3  ro 

1 

o 

1 

o o 

1 

o 

1 t 

0 0 

1 1 

0 0 

1 

0 

0 

1 

CO 

0 0 

0 

000 

o 

(O 

CO 

CO  to 

to 

>0 

cn 

ro 

XU 

XU 

XU 

0 ro  to 

-< 

tn 

ro  o 

o 

— # CD 

XU 

0 

ro 

CD 

—* 

cn 

ai 

CO  to  '*0 

X 

<_n 

ro 

cn 

Oy  — ' 

cy 

to 

CD 

cn 

CO 

CO 

ay 

CT  OJ  0 

v: 

f\j 

ro 

CD 

(0 

•*0 

CO 

ro 

CD 

XU 

ro 

jy  ro  CO 

o 

0 

' 

1 

1 

1 

1 

' 

1 1 

D< 

o 

o 

o o o 

0 

0 

CD  0 ro 

0 

0 

0 

0 CD  0 

33 

0 — ' 

0 

0 

0 

*-o 

CD 

to 

x» 

r>o  0 — ^ 

-< 

<_n 

ro 

to  U3 

iO 

to 

0 

rn 

XU 

ro 

CD 

XU 

CO 

t->  to 

43 

cn  00 

CO 

CO 

cn 

ro 

.XU 

cn 

ro 

'o  0 <n 

rsi 

0 

CO  o 

CO 

C3 

•'“* 

cn 

<7! 

VO 

tn 

ro 

to 

•^4 

ro  — • VO 

X 

-H 

1 

1 

1 

» 1 

1 

1 

1 

1 

1 

1 1 1 

r“ 

• 

• 

• 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• • • 

0 

o 

o o 

0 0 

0 0 

0 0 

0 0 

.mm* 

000 

0 

<-o  ay 

oy 

cn  cn 

£u 

0 

lu 

XU 

Xu 

-c*  0 

Xu  — * ro 

•< 

cry 

ay 

03 

CO  ro 

cn 

CO 

-t* 

<33 

ro 

CD  CD  to 

VO  CO  0 

M 

Ln 

ro  o 

cn 

ro 

CD 

cn 

-0 

ro 

CO 

00 

cn 

VO  ro  ro 

X 

o 

CJ>  oy 

rvj 

ro 

to 

0 

CO 

CO 

cn 

0 CO  0 

o 

rr 

Q. 

3 

vO 


3" 


3 


TD 

O 

o 

o 

l/> 

trv 


J 


t 


Table  3.  Change  in  engineering  shearing  strains  in  the  6 regions  and  17  tetrahedrons  of  the 
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all  of  the  tetrahedrons  which  share  the  node  under  consideration. 
For  example,  node  1 is  shared  by  the  tetrahedrons  (1),  (2),  (3), 

(10)  and  (11)  as  indicated  in  Figure  13.  Thus,  as  far  as  at 

node  1 during  the  first  deflation  step  (the  column  of  100%  80% 

TL^C  in  Table  2)  is  concerned,  the  values  -.03718,  -.13338, 

-.04777,  -.04535,  and  -.04651  for  the  five  sharing  tetrahedrons 
can  be  combined  to  give  an  average  value  equal  to  -.06204.  This 
process  of  averaging  has  been  applied  for  all  nodes  and  the  re- 
sults are  presented  in  Table  4 (see  page  3 4)  . For  interpretation 
of  the  strain  results,  the  data  in  Tables  2 and  3 are  to  be  used 
for  the  strain  changes  at  the  centroids  of  the  tetrahedrons,  and 
Table  4 for  those  changes  occurring  at  the  nodes. 

When  changes  are  plotted  (Figure  14) , it  is  observed  that  the 
lower  layers  of  both  regions  A and  B respond  with  smaller  than 
the  other  layers  at  the  beginning  stage  of  the  volume  decrease  but 
tend  to  follow  the  general  trend  of  the  other  layers  at  later 
stages  of  the  deflation  process.  It  is  also  noteworthy  that  the 
strain  variations  can  also  be  studied  around  the  surface  of  the 
lobe  by  traversing  along  the  lines  1-5-13-18,  2-6-14-19,  2-9-12-20, 
and  3-4-10-16  which  are  the  four  partitioning  lines  used  in  the 
formation  of  the  tetrahedrons.  The  information  needed  for  such 
study  can  be  extracted  from  Figure  14.  For  example,  the  strain 
data  for  the  line  1-5-13-18  are  those  illustrated  in  the  second 
column  of  each  of  Figures  14(a)  through  14(h);  the  node-1  column 
of  Figure  14(e)  may  also  be  considered  as  second  column  with  the 
first  and  fourth  columns  unavailable. 


Changes  of  e are  found  similar  to  those  of  e . 

y ^ 


The  apex-to-base 


variations  of  strains,  e , are  generally  greater  in  magnitude 

z 

than  and  e^.  And  this  disparity  between  the  radial  strains, 

and  e^,  and  the  axial  strain,  is  more  pronounced  at  the 

basal  (lower)  region  than  at  the  apical  (upper)  region  of  the 

lobe.  Figure  15  shows  the  relatively  uniform  region-by-region 

and  layer-by-layer  e_  changes  in  the  excised  lobe. 

z 
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Table  4.  Changes  in  the  average  cnninccrinfi  strains  at  tfie  nodes 


(metallic  markers)  during  tlic  deflation  process. 


Node 

No. 

% TL  C 

0 

^7 

Y 

' XV 

100^80 

-.06204 

-.07037 

-.08901 

-.00448 

-.01483 

.01571 

1 

10CH60 

-.12282 

-.1  5052 

-.16439 

.00847 

-.01511 

-.00745 

lOCMO 

-.18651 

-.21484 

-.25693 

-.00815 

-.00019 

-.04191 

1 OO-^^BO 

- .04656 

-.06676 

-.08548 

.02334 

-.00759 

-.00619 

2 

100-60 

-.12561 

-.15482 

-.16387 

.01367 

-.02434 

-.00775 

100-40 

-.17089 

-.20358 

-.25525 

-.03326 

.01479 

- .06192 

100-80 

-.09058 

- .09199 

-.10280 

-.02365 

.00182 

.06209 

3 

100-60 

-.13179 

-.16522 

-.18089 

.00501 

.0178“ 

.0415? 

100-40 

-.15670 

-.24341 

-.26306 

.01359 

.01036 

-.03460 

100-80 

-.07062 

-.08350 

-.10042 

.01917 

-.01119 

.04398 

4 

lOO-'^O 

-.12160 

-.15524 

-.19959 

.02280 

-.01000 

.02667 

100-40 

-.19387 

-.24386 

-.29262 

.03579 

-.01781 

.00347 

100-80 

-.03917 

-.06657 

-.07694 

.01860 

-.00681 

-.01907 

5 

100-60 

-.10393 

-.14801 

-.17739 

.01794 

-.01788 

-.03476 

100-40' 

-.21107 

-.21742 

-.28229 

.01152 

-.01137 

-.02341 

100-80 

-.04293 

-.04182 

-.07044 

-.01401 

-.01346 

-.01976 

6 

100-60 

-.11578 

-.12755 

-.16419 

-.01654 

-.03920 

-.04146 

100-40 

-.21352 

-.17560 

-.26016 

-.04396 

-.00432 

-.05893 

1 00-80 

-.04249 

-.0638? 

-.07743 

.00035 

-.01368 

.01866 

9 

100-60 

-.10748 

-.14241 

-.16466 

.,00660 

-.02013 

.00612 

100-40 

-.19117 

-.21369 

-.25938 

-.00016 

-.01465 

-.00468 

100-80 

- .00324 

- .06068 

-.07312 

.00097 

.00299 

.01162 

10 

100-60' 

-.05697 

-.12429 

-.17419 

-.02144 

.01626 

.00242 

100-40 

-.14319 

-.19733 

-.29147 

-.04041 

.00793 

.00230 

100-80 

-.00227 

-.05746 

-.05155 

-.00053 

-.03508 

.06247 

12 

100-60 

-.06773 

-.13170 

-.13370 

-.00581 

-.03595 

.07905 

100-40 

-.15982 

-.21602 

-.25619 

-.00366 

- . 06096 

.09696 

100-80 

-.01204 

- .04504 

-.05714 

-.01211 

-.02175 

.03306 

13 

100*60 

-.07022 

-.1  0479 

-.14881 

-.02785 

-.02428 

. 02871 

100-40 

-.15270 

-.17472 

-.25750 

-.03405 

-.02903 

.00221 

100*-80 

-.00792 

-.0355? 

-.05214 

- .03128 

-.03988 

.0^1569 

14 

1 00-60 

-.07496 

-.09821 

-.14349 

-.03203 

-.05944 

.05463 

100-40 

-.16793 

- . 1 81 38 

-.25175 

-.02003 

-.06000 

.03397 

100-80 

.00778 

-.01325 

-.06827 

-.05441 

-.00897 

.04078 

16 

100-60 

-.01720 

-.02116 

-.15901 

-.14629 

-.03964 

.00528 

100-40 

-.06581 

-.06267 

-.27367 

-.22031 

-.05640 

-.07343 

100-80 

-.04202 

-.03713 

-.0461“ 

.01458 

-.05377 

.0“139 

18 

1 00>-60 

-.09665 

-.07601 

-.13829 

-.00389 

-.07973 

.09638 

100-40 

-.14340 

-.13446 

-.22418 

-.03723 

-.08576 

-.01054 

100-80 

-.01 500 

-.03907 

-.05130 

-.00559 

-.03884 

.05704 

19 

100-60 

-.07648 

-.09149 

-.13814 

-.00402 

-.07300 

.06547 

100-40 

-.1  3767 

-.16512 

-.22511 

.00338 

-.09648 

-.02638 

100-80 

-.00276 

-.03679 

-.05417 

-.02111 

-.02683 

.04251 

20 

100-60 

-.05352 

-.08223 

-.1504? 

-.04849 

-.02745 

.04253 

100-40 

-.12059 

-.14779 

-.26996 

-.07245 

-.03196 

-.01398 
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As  shown  in  Table  4 and  Figure  16,  the  shearing  strains  are  small 
as  compared  with  the  normal  strains  except  in  the  vicinity  of 
the  node  16.  The  reason  for  such  behavior  of  the  lobe  in  this 
vicinity  being  different  from  other  localities  could  be  the  near- 
ness of  this  node  16  to  the  cannula  which  was  inserted  into  the 
main  bronchus  of  the  lobe  and  held  it  securely  during  the  de- 
flation experiment. 

The  avei ige  strains  in  layers  (I),  (II),  and  (III),  and  in  the 
regions  A and  B can  be  obtained  by  combining  the  strain  data  for 
the  tetrahedral  elements  in  the  layers  and  the  regions,  respectively. 
The  results  are  listed  in  Table  6.  Generally  spealcing,  region  A 
exhibits  greater  normal  strains  but  smaller  shearing  strains 
than  region  B,  and  the  basal  layer  (III)  exhibits  smaller  normal 
strains  but  greater  shearing  strains  than  the  apical  layer  (I) 
in  both  regions  A and  B.  Since  region  A and  basal  layers  contain 
relatively  more  bronchi  than  region  B and  apical  layers,  such 
disparities  in  the  strains  during  the  inf lation -deflation  process 
might  be  expected. 

Displacement  functions  other  than  those  given  in  Equations  (1)  and 
large  strain  equations,  instead  of  small  strain  equations  (Equations 
(4)),  are  being  studied  so  that  the  developed  method  can  be  expanded 
and  applied  not  only  for  the  excised  lung  lobes  but  also  for  the 
intact  worlcing  lungs,  and  furthermore  for  the  evaluation  of  the 
rotation  and  tra.islat  on  of  the  intact  wor)cing  lungs  during  various 
respiratory  maneuvers  a i thoraco-abdominal  configurations. 
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Overall  Program  Summary  and  Significance.  The  development  in  our 
laboratory  of  the  parenchymal  marker  technique  utilized  in  Parts  I 
and  II  as  well  as  the  ability  to  reconstruct  the  three-dimensional 
shapes  and  dimensions  of  the  heart,  lungs,  chestwall  and  diaphragm, 
in  the  intact  thorax  (detailed  in  last  year's  progress  report) 
presents  a unique  opportunity  to  study  dynamic  regional  lung 
mechanics  and  lung-chest  wall-diaphragm  interaction  under  various 
conditions  of  lung  volumes  and  body  position  and  thereby  provides 
the  methodology  to  study  the  effects  on  regional  lung  function  of 
alterations  in  these  parameters  induced  by  changes  in  the  gravita- 
tional-inertial force  environment.  The  results  described  in  Parts 
I and  II  represent  significant  new  information  regarding  regional 
mechanical  properties  of  the  lung  and  clearly  indicate  the  need  to 
study  the  intrinsic  elastic  behavior  of  the  lung  and  its  relation- 
ship to  the  chest  wall  and  diaphragm  in  the  intact  thorax  with 
techniques  possessing  sufficient  temporal  and  spatial  resolution 
to  permit  quantitative  determination  of  the  dynamic  changes  in 
spatial  (three-dimensional)  lung  parenchymal  strains,  regional 
lung  volumes  and  lung-chest  wall  geometries. 

Mechanical  studies  of  the  intact  lung  have,  heretofore,  not  been 
possible.  These  data,  which  cannot  be  obtained  by  other  methods, 
provide  a quantitative  measure  of  the  mechanical  properties  of 
the  lung  under  normal  physiologic  conditions  (i.e.,  an  intact 
thorax  and  normal  intrathoracic  pressures)  at  IG  and  are  required 
before  any  similar  studies  can  be  performed  under  conditions  of 
increased  or  decreased  gravitational-inertial  force  environments 
encountered  in  aerospace  flight. 

Because  of  the  very  large  differences  in  specific  gravity  of 
the  air  in  the  alveoli,  and  the  lung  parenchyma  and  its  mobile 
blood  content  plus  the  anatomically  fragile  nature  of  the  lungs, 
the  respiratory  system  is  the  most  susceptible  of  all  bodily 
organs  to  malfunction  and  actual  structural  damage  during  exposure 
to  changes  in  the  direction  and  magnitude  of  the  gravitational- 
inertial  force  environment. 
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Consequently , increased  knowledge  concerning  the  mechanical 
characteristics  and  dynamic  three-dimensional  geometry  of  the  lung 
is  a basic  requirement  for  understanding  the  factors  which  deter- 
mine man's  reactions  and  tolerance  to  short  and  long  term  changes 
in  the  force  environment  such  as  encountered  in  conventional  and 
high-G  aircraft  and  in  manned  space  flight. 


Project  Director 


Peter  A.  Chevalier,  Ph.D. 
Principal  Investigator 


November  23,  1977 
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